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Active magnetic bearings can be implemented with frequency band-reject filtering 
that decreases the bearing stiffness and damping at a small bandwidth around a chosen 
frequency. The control scheme has been used for reducing a rotor dynamic force, such as 
an imbalance force, transmitted to the bearing stator. This study reveals that the 
scheme creates additional system vibration modes at the same frequency. It also shows 
that the amount of force reduction is limited by the stability requirement of these 
modes . 


INTRODUCTION 

The attractive-type active magnetic bearings (AMBs) usually have four quadrants of 
electromagnets (ref. 1). A pair of opposite quadrants independently control the jour- 
nal motion in one direction by a Proportional-Integral-Derivative (PID) controller 
(ref. 2). Steady-state (bias) currents are induced in the quadrants so that the total 
control current in each quadrant never changes polarity. This provides a base for line- 
ar feedback control and greatly simplifies the control circuitry (ref. 3). 

There appears to be a growing interest in applying the AMBs not only for their 
advantageous basic bearing function, but also for their potential to be the rotor force 
isolator (ref. 4). An AMB can be made extra "soft" at narrow frequency bands, and the 
rotor forces to ground at these frequencies can be dramatically reduced at the sacrifice 
of large rotor runouts. For example, the well-publicized AMB control feature, 
'’auto-balancing, 11 was designed to reduce imbalance force to ground with the rotor 
rotating about its inertia axis (ref. 2). However, the creation of these stiffness 
"valleys" also creates instability problems for the rotor-AMB system. Herein this type 
of instability problem is addressed analytically while an experimental project is in 
progress. 

NOMENCLATURE 


AMB active magnetic bearing 

a phase-lead network zero parameter 

B frequency bandwidth of band-pass filter 

b phase-lead network pole parameter 

C radial air gap 

c d proportional feedback gain 
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C v phase-lead feedback gain 

C e integral feedback gain 

Cp band-pass filter gain 

Cp 0 band-pass filter stability threshold gain 

F regulating magnetic force 

^ 1*^3 magnetic force per pole due to bias current in quadrant 1 or 3 

FODE first-order differential equation 

f band-pass factor 

G power amplifier gain 

Ij.,13 bias current in quadrant 1 or 3 

i regulating current 

j /=T 

K AMB stiffness 

current stiffness 

K m magnetic (negative) stiffness 

LHS left hand side 

M mass supported by AMB 

PID proportional-integral-derivative 

Q integrator output 

R nondimensional frequency parameter defined at a stiffness valley 

RHS right hand side 

S Laplace variable 

t time 

U real stiffness at U) c after filter implementation 

u real stiffness at u> c before filter implementation 

u r RHS stiffness 

ufl, LHS stiffness 

V imaginary stiffness at U) c after filter implementation 

v imaginary stiffness at U) c before filter implementation 

Y AMB journal displacement in Y-direction 

Yp band-pass filter output 

Z phase-lead network output 

ACrf net proportional feedback gain contributing to real stiffness 

U) exciting frequency 

0) c band-pass filter center frequency 

U)^ LHS slope frequency 

U) q integrator cut-off frequency 


342 



U)^ RHS slope frequency 

dot d/dt 

Au> 0) -U) 

c 

0 magnetic pole angle 


AMB CONTROL WITH FREQUENCY BAND FILTERING 

Assuming the bias currents are much larger than the regulating currents and the air 
gap is much larger than the journals normal excursion, the two perpendicular axes of a 
radial AMB can be controlled independently (ref. 3). A control scheme for each axis is 
presented in figure 1. There are four parallel loops processing a journal displacement 
measurement. The top three loops, including a phase-lead network (b > a), form a 
conventional PID controller. The fourth loop comprises a typical second-order 
band-pass filter with the center frequency at bandwidth B, and gain Cp (ref. 5). 
The filter can have a fixed center frequency or it may vary the center frequency with 
rotor speed. The latter is called a tracking filter and its implementation was 
explained by SKF (ref. 6 ). The filter output is subtracted from the displacement meas- 
urement that is fed into the phase-lead network. Also, the output multiplied by a gain 
ACrf is added to the basic PID signals that control the power amplifier. The gain ACj is 
defined below. 

Ac d = C d - K m /(KiG) (1) 


where 


Cd = proportional loop gain, volt/m 

G = amplifier gain, A/volt 

K m = magnetic stiffness 

= 2 (Fj/C + F3/C)cos 0, lb/m 

= current stiffness 
= 2 (Fi/Ii + F 3 /I 3 ), lb/ amp 

C = radial air gap, m 

Il,l 3 “ bias current of quadrants 1 or 3, respectively, A 

Fl* F 3 = magnetic forces per pole due to bias currents, N 

The parameter K m is the "negative spring" effect of the AMB magnetic field. The 
part of the proportional gain to overcome this effect is K m /(K£G). Therefore, ACj is 
the net proportional gain contributing to the AMB stiffness. 

Assuming the power amplifiers are current sources in the bandwidth of interest (G = 
constant), the AMB regulating force across the quadrants 1 and 3 is 

F = Ki i + K m Y (2) 

From figure 1, the regulating current is 

i * G [ - C d Y - C v Z - C e Q + Ac d C p Y p ] (3) 
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where 


z = [(S + a)/(S + b)] (Y - CpYp) 

Q = [1/(S + W Q )] Y 

Y d = [BS/(S 2 + BS + U) 2 )] Y 


(4) 

(5) 

( 6 ) 


Incorporating equations (4), (5), and (6) into equation (3), which in turn is incorpo- 
rated into equation (2), the AMB reaction transfer function can be expressed by equation 
(7). 


-F/Y = KiG[(Ac d + C V (S + a)/(S + b))f + C e /(S + U) o )] 


(7) 


where 


f = 1 - C d BS/(S 2 + BS + 0) 2 ) 

r C 

The filter gain Cp ranges from 0 to 1. For AMB without frequency band rejection, 
i.e., Cp - 0, the complex stiffness as a function of exciting frequency is 

-F/Y = [K£G(Ac d + C v (ab + U) 2 )/(b 2 + oj 2 )) + KiGC e u) o /(u) 2 + W Q 2 )] 

+ j 0)[KiGC v (b - a)/(b 2 + u> 2 ) - KiGC e /(to 2 + U) q 2 )] (8) 

The real part of equation (8) is the AMB stiffness and the imaginary part is the 
AMB damping. The second terms of both parts are the main contributors to the AMB "stat- 
ic stiffness" (for small to). The stiffness and damping of a typical AMB are presented 
in figure 2. It may appear unusual that the bearing damping can be negative, but as 
long as no rotor natural vibration mode exists in the frequency range with negative 
damping, there should be no dynamic problem. 

Figure 3 shows the magnitude of the complex stiffnesses for the same AMB with a 
frequency band rejection. The filter creates a stiffness valley with a depth propor- 
tional to Cp. WithCp approaching one, adynamic force exerted on the rotor at the center 
frequency can be mostly balanced by the rotor inertia force. Only a small part will be 
resisted by the AMB and thus transmitted to the bearing stator or ground. In the 
following section, a potential stability problem of creating such a stiffness valley 
will be discussed. 


AMB STABILITY AT FILTER CENTER FREQUENCY 

Figure 4 is a zoomed-in view of the stiffness valley of figure 3. The local 
complex stiffness decreases and increases sharply but continuously around a center 
frequency to c in a small bandwidth B. 

Let 


M = rotor mass associated with the AMB 

u, v = real and imaginary parts of complex stiffness, respectively, 
at before the valley was created 
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and 


/uTM > U) 

c 

which implies a natural frequency exists above the filter center frequency. 

Imagine applying a dynamic force to the mass at a frequency co that is slowly 
increasing across the valley. There will be (if the valley is deep enough) an exciting 
frequency, U)^, associated with a LHS slope stiffness, u£, such that, 

ATjjM - U)£ 

Therefore, is a resonance frequency. Similarly, there is a RHS slope resonance 
frequency: 

/77m - a) 

r r 

These discussions are best illustrated by an example shown in figure 5. Note that 
these two resonance modes are not due to the filter circuitry alone. They are also 
related to the mass which the AMB sees. It is of less concern how well the modes are 
damped. Presumably, if the exciting force frequency, such as the rotor speed, is not 
drifting away from the filter center frequency, these modes can only be excited by 
impact type loads. Since the filters are usually implemented with narrow bandwidth, 
they will not be excited easily as long as they are reasonably damped and no persistant 
impact load exists. It is a major concern, however, that these modes may not be stable, 
i.e., associated with positive growth factor or negative damping. This can happen when 
the gain Cp is made large or approaching one. For example, the LHS slope mode at 59.4 Hz 
in figure 5 is unstable (growth factor = 120) with Cp = 1. It is therefore important in 
implementing this type of filter to know a threshold stable gain Cp, which is determined 
in the following analysis. 


Let 

and 


Au) * a) - a) 

c 

R = (B/2)/Au> 


The band-reject factor in equation (7) is 


f = 1 - C p (jBu>)/[(u) c 2 - U) z ) + jBw] 


“ [1 + (1 * C p )R 2 - jC p R]/(l + R 2 ) 

(9) 

The complex stiffness around U) c is approximately 


-F/Y « ( u + jv) f = U + jV 

(10) 

where 


U = u [1 - C p R Z /(l + R 2 )] + V CpR/(l + R 2 ) 

(11) 

V = v [1 - C p R 2 / ( 1 + R 2 )] - u CpR/(l + R 2 ) 

(12) 
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Both U and V are functions of R, which in turn is a function of the exciting 
frequency near 0) . According to equation (12), the RHS mode (R < 0) is always more 
damped than the LHS mode. Note that at the center frequency (R = °°), equations (11) and 
(12) become 


U = u ( 1-Cp) 

and 

V = v (1-Cp) 

Also note that for effective force isolation, the frequency /[u(l - Cn)]/M should 
be one-third of w c or less. To determine the threshold gain (Cp Q ) for stability, it may 
not be overly conservative to require 

V > 0 

which implies by equation (12) 

C po < (1 + R 2 )/[(u/v)R + R 2 ] (13) 

In the range of 0 to 1, the minimum value of Cp Q occurs at the LHS slope where 

R = v/u + /v 2 /u 2 + 1 (14) 

Incorporating equation (14) into equation (13), the relationship between Cp Q and v/u is 
plotted in figure 6. For normal AMB applications with v/u < 1, the gain value of Cp 
should be less than 0.83 according to this plot. 

STABILITY OF ROTOR-AMB SYSTEM USING FREQUENCY BAND FILTERING 

When two or more radial AMBs are supporting a rotor, the mass that each AMB sees is 
different at different critical modes. The location of the filter center frequency 
relative to these critical frequencies has a definite influence on the stability prob- 
lem mentioned above. Since the influence is not straight forward, it would be appropri- 
ate to investigate the stability problem in a rotor-AMB dynamic system as follows. 

In a conventional rotordynamics approach, the rotor is modeled as sections of 
circular beams using a finite element method. Concentrated masses and inertias are 
assigned at the nodes of the beam elements for any attachments to the rotor. Gyroscopic 
effect is included. For simplicity, circular orbits can be assumed and are adequate for 
most AMB amplications. For each radial AMB there are two independently controlled axes. 
For each axis, there is a set of first-order differential equations (FODEs) represent- 
ing the AMB dynamics. For example, the control scheme of figure 1 can be represented by 
four FODEs according to equations (4), (5), and (6), which include the frequency band 
filtering. The rotor-bearing coupling terms exist in equations (2) and (3). 

Combining the dynamic equations of the rotor and AMBs, an electromechanical system 
model can be formulated for eigenvalue evaluation. The formulation procedure is 
straight forward and will not be presented here. An example eigenvalue analysis for a 
simple rotor supported by two identical AMBs (fig. 7) is used to demonstrate this system 
approach. The key AMB control parameters are identical for all four controlled axes 
(table I). The stiffness and damping of this AMB as functions of frequency have been 
plotted in figure 2. 
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Without the filtering, the system has natural modes at 3000 cpm, 10,000 cpm and 
35,000 cpm. The rotor-AMB system has been analyzed for the filter centered at 30 Hz and 
60 Hz separately. The additional LHS mode, which is less stable than the RHS one, is 
presented in table II for different values of Cp. The predicted values of Cp 0 using 
equations (13) and (14), as noted on table II, are consistently conservative. 

CONCLUSIONS 

A frequency band-reject scheme for reducing dynamic force to stator at a selected 
frequency may create AMB instability problems. 

From the generic bearing point of view, this study has shown the following: 

1. Depending on the mass supported by the AMB, there can be a natural vibration 
mode corresponding to the stiffness somewhere on each slope of the stiffness 
valley. Thus the frequencies of the induced modes are close to the filter 
center frequency. 

2. There is a limit how deep the stiffness valley can be made without causing 
these modes to be unstable. The limit is related to the local damping-to- 
stiffness ratio before the valley is created. A conservative limit in terms of 
the ratio has been established. 

From the rotor-AMB system point of view, the AMB mass varies with the rotor crit- 
ical mode shapes. A rigorous approach to determine the stability is to find the eigen- 
values of the electromechanical system. 
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TABLE I - RADIAL AMB PARAMETERS 


Ki 53.5 N/A 

5 

K m 2.27x10 N/m 

a 163.3 rad/sec 

b 978.8 rad/sec 

to 3.14 rad/sec 

o 

Cd 5921 volt/m 

C v 5626 volt/m 

C e 31500 volt/m-sec 

Cp 0 to 1 

to 188.5 rad/sec or 377.0 rad/sec 

c 

B 7.54 rad/sec or 15.08 rad/sec 


TABLE II - ADDITIONAL LHS MODE 



Speed = 

O) = 
c 

v/u = 
Cpo = 

1800 rpm 
30 Hz; 

0.25; 

0.4 



Speed = 

O) = 
c 

v/u = 
Cpo = 

3600 rpm 
60 Hz; 
0.46; 

0.6 

c p 

Frequency log, 

(cpm) decrement 

Frequency log. 

(cpm) decrement 

0.0 

1800 

0.126 

3600 

0.126 

0.4 

1798 

0.045 

3597 

0.064 

0.5 

1798 

0.025 

3597 

0.049 

0.6 

1797 

0.006 * 

3596 

0.034 

0.8 

1796 

-0.034 * 

3595 

0.003 * 

1.0 

1794 

-0.072 * 

3593 

-0.027 * 


* marginal or unstable modes 
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AN ACTIVE MAGNETIC BEARING CONTROL SCHEME 
WITH FREQUENCY BAND-REJECT FILTERING 



Bias I3 
Adjustment 


881081-1 


Figure 1 


STIFFNESS AND DAMPING OF A TYPICAL ACTIVE MAGNETIC BEARING 



Figure 2 


349 
















Natural Frequency (Hz) 


NATURAL FREQUENCY VS. EXCITING FREQUENCY 


Local AMB 
Stiffness and 
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881491 

Figure 5 


C po VS v/u 



v/u 

881494 

Figure 6 
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ROTOR MODEL 
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Figure 7 
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